Dynamical Mean Field Theory of Temperature and Field Dependent Band Shifts in 
Magnetically Coupled Semimetals: Application to EuB 6 
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A model for semimetals such as EuBe, in which band overlaps are controlled by magnetic order, 
is presented and is solved in the dynamical mean field approximation. First order phase boundaries 
are computed by evaluating free energies of different states. The phase diagram is determined. A 
specific and physically reasonable choice of parameters is found to approximately reproduce the 
available data on EuBq. For this material, predictions are made for the location of a metamagnetic 
transition and its associated endpoint, and a change in the order of the magnetic transition. 
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I. INTRODUCTION 



II. MODEL HAMILTONIAN 
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EuBq is a magnetic semiconductor in which magnetic 
order apparently controls the bandgap. The ferromag- 
netic transition at T c ~ 12K is signaled by a sharp 
increase in plasma frequency and a sharp drop in 
resistivity Q. In this paper we formulate a model 
which captures the physics of EuB e , and solve it in the 
dynamical mean field approximation. The important 
technical step is the accurate determination of first order 
phase boundaries via the construction of the free ener- 
gies of competing phases. We present a general phase 
diagram, which includes a change in the order of the 
transition from first to second. Remarkably, the change 
occurs not by the vanishing of a fourth order coefficient 
but simply by an exchange of stabilities. We show that 
for a choice of parameters consistent with what is cal- 
culated for EuBq 0, the various experimental data are 
semiquantitatively reproduced. For these parameters, we 
predict that EuB 6 should exhibit a metamagnetic phase 
transition, and we estimate the location of the transi- 
tion line and its endpoints. The change in order of the 
phase transition might be accessible in pressure exper- 
iments. The analysis reported here may be thought of 
as a more precise solution of a model proposed by 
Our solution agrees with theirs in essential aspects but 
provides a sharper picture of physics and includes some 
new features. 



We also note that the band structure of the material 
is presently the subject of controversy with photoemis- 
sion experiments indicating a large gap (implying that 
the material is intrinsically insulating so the metallic be- 
havior is due to defects which dope the system) while 
some band calculations Q [(| and quantum oscillation 
measurements 8J indicating a small negative gap causing 
the semimatal behavior. We comment below on the im- 
plications of our work for these controversies. 



The electronic structure of EuBq has been calculated 
[j| 0. Experiment, quantum chemical intuition and 
band calculations all agree that the Eu /-shell is half 
filled and, as the scale relevant for electronic behavior, 
electronically inert. Crystal field effects are negligible 
and to a good approximation, each Eu may be regarded 
as carrying a S — 7/2 "core spin" . Band theory calcu- 
lations jfj reveal two near Fermi surface bands: a nearly 
empty band with a minimum at X-point, derived mainly 
from Eu — d, B — p orbitals, and a nearly full band with 
a maximum at the X-point, derived mainly from Eu — /, 
B — p ones. The hybridization between these bands is 
negligible because they arise from different symmetry or- 
bitals. 

Because we shall be interested in low energies and weak 
couplings, we expand the bands near the X-point so the 
minimal model describing EuB§ becomes: 
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Here S represents the S = 7/2 Eu core spin, c\ and ci 
represent upper (conduction) and lower (valence) bands 
and A controls the band overlap. We take \S\ = 1, ab- 
sorbing the actual magnitude into the coupling J. Be- 
cause S originates from a filled (spin polarized) Eu shell, 
and the latest band structure calculation ;0j shows that 
two bands arises from different Eu orbitals (/ and d), we 
may therefore take J to be diagonal in orbital indices. 
The high spin and low carrier density means the Kondo 
effect is irrelevant, so the sign of J is arbitrary. Kunes 
and Pickett's calculation 6] implies J 11 < (Kondo cou- 
pling) and J 22 > (anti-Kondo coupling) and indicates 
that the splitting of the conduction and valence bands 
are roughly the same, we choose J = J 22 = — J 11 . Our 
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FIG. 1: Expanded view of ferromagnetic phase band structure 
of EuBe with momentum k measured from X [n, 0, 0], in units 
of inverse lattice constant a — 4.2A Curved (flat) dashed line 
represents valence (defect) band and heavy line conduction 
band. 

results turn out to depend mainly on |J 22 | + |</ n |. The 
chemical potential fi is determined by the neutrality of 
the system, i.e. number of particles equals number of 
holes. FigQJ shows the simplified band structure implied 
by eqnJTJ. 

The model described in eqn([Q) can be easily modified 
to match the bandstructure implied by the photoemission 
experiment as follows: we interpret the lower band (02) 
as arising from defect states which in this picture must 
exist so that carriers are donated to the conduction band, 
and thus take m-i to be infinite, and assure there is no 
coupling between C2 band and the core spin S (J 22 = 0). 
At these parameters (ma — > 00, J 22 = 0), the C2 band 
serves as a particle reservoir and the chemical potential 
is fixed at 0. This choice in fact does not change the 
structure of the phase diagram or the qualitative features 
of our results. 

The T — phase diagram is straightforward. The 
ground state is ferromagnetic. If A + \ J 11 \ + \ J 22 \ < 0, 
the majority spin of band 1 crosses the minority spin 
of band 2, which results in metallic behavior. If not, 
the ground state is insulating. At T > T c , the random 
spin orientations mean that the band shift are much less 
(especially at small J). Roughly, A < implies metallic 
behavior while A > implies insulating state. 



III. BASIC RESULTS 
A. Dynamical Mean Field and Free Energy 

To study eqnl^Q), we use the single site dynamical mean 
field method 9] . This amounts to assuming that the elec- 
tron self energy is momentum independent: £ 1-0 .(p, ui) — > 
Ei i(T (w). The self energy is determined from the so- 



lution of an auxiliary impurity model with free energy 
^imp — ~Tlog Zi mp , and 

Z imp = jDlc+c]D[S c ]ex P S eff 

with [13 

1 f 
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The free energy of the lattice model is 
= Qimp - ry^log G a (iw n ) 

J deD ( e ) x l og[iLu n +n- S CT (iw n ) - e] (3) 

n.cr 

and the auxiliary function a is fixed by the requirement 
^ = 0. It is important to perform the exact quantum 
trace over the spin degrees of freedom rather than em- 
ploying the classical spin approximation often made 01 
in order to obtain physically reasonable estimates for the 
entropy. 

The free energy so determined is a functional of the 
applied field h and the system magnetization is given by 
m = — |P. It is sometimes useful to perform a Legendre 
transformation to obtain f2(m, T). We have not found 
an efficient method for performing this transformation: 
to obtain n(m,T), we solve the dynamical mean field 
equations numerically for a range of h, and then con- 
struct f2(m, T) = f2(/i, T) + h m explicitly. We note how- 
ever that in some parameter ranges, the DMFT equa- 
tions have two stable (symmetry un-related) solutions at 
h = 0, corresponding to states with m — and m ^ 0. 
The two states so determined are extrema of f2(m, T) and 
their free energies are the extremal values of Q(m, T), so 
first order transition points may be located without con- 
structing f2(m, T) explicitly. 

B. Approximations and DMFT Solutions 

We have solved the dynamical mean field equations 
corresponding to extremizing eqnQ with respect to 
S(w„). For simplicity we adopted a semicircular den- 
sity of states Die) — ^^j/* with t = < '^ 2lT X s chosen to 
match the band theory band-edge density of states. We 
also assumed J/t « 1 (because this is the limit relevent 
for the materials of interest) and T/J<< 1 (we will see 
below that the calculated transition temperature is much 
less than J) . These approximations simplify calculations 
considerably, in particular, in this limit we may retain 
only the z component of the core spin, simplifying the 
quantum trace. 
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To simplify the expression, we define 
A = 2yJ 



(4) 



where y is a dimensionless parameter measuring the 
bandgap. We also define the following / functions: 

f 15 ( x ) = e(x)x 1 - 5 
f 5 (x) = e(x)x 2 - 5 



where Q(x) is the step function. 
In the small J limit, 

S T,(1) = T mJ 



(5) 



the free energy Af2 and core spin magnetization m are 
given by 
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m = ~ sinh( y ) - \- sinh(~ ) = S 7/2 [a (to, /3)] (7) 



where 



a = — (3 J x 

6TT 



■) 1J (/ 15 (M - (V - m)) - / 15 ( M - (to + y))) + 
(£) (/ 15 (m - y - /x) - / 15 (-(m + y) - fij) (8) 



The chemical potential /i is determined by charge neu- 
trality, i.e. number of electron equals number of holes. 

J_ (/15 (jU _ (m _ y)) + / lB (/i _ (m + 

= 7T5 (/ 15 (™-2/-M) + / 15 (-(^ + ?y)-M)) (9) 

Note that is a function of m and y, not an independent 
variable. 

Based on the above expressions, the entropy S is 
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the specific heat Cy is 

d AE d(An + TS) 



dT 
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FIG. 2: The phase diagram. In the range — 0.1J < A < 
— 0.55J, the system undergoes both 1st and 2nd order phase 
transitions. The temperature is in unit ^rii(-p) 2 ' 5 . To = 
17.43 if if the J = Q.lGeV and t\ = lA6eV implied in 
Kunes/Pickett's bandstructure calculation. 



can be computed from eqnJ7|). 
Finally, the density of conducting band electrons An 

is 

An = ) 3 / 2 (/ 15 (u(m, y) - (y - m)) + 

/ 15 (u(m, y) - (m + y)) - 2 / 15 ( M (0, y) - y)) (12) 

The plasma frequency w p is defined as = 47re 2 and 
calculated by 



<4(T-)-wl{T+) 



47re 



, An 



(13) 



We emphasize that all expressions above should be 
evaluated at the m where eqn(J2J| is satisfied. 

The calculated phase diagram is shown in Fig©. We 
see, in accord with the simple considerations of the pre- 
vious section that for A/ J > 1, the material is always 
insulating and T c is negligible. For A sufficiently nega- 
tive, the transition becomes second order. 

Fig© shows the the free energy as a function of m at 
several temperatures and at different bandgaps. Panel 
(a) shows f2(m) for several T at A/J = 0, in the first 
order region of the phase diagram. Panel b shows the 
unusual behavior in the vicinity of the multicritical point 
where the transition becomes second order. One naively 
expects the minimum at higher m shifts downwards, 
eventually merging with m = minimum (in other word, 
that the sign of the m 4 term in the Landau expansion of 
the free energy changes) . Panel (b) shows that the actual 
situation is more subtle: the higher m minimum contin- 
ues to exist; however its energy is increased so the second 
order phase transition happens, and then at a lower tem- 
perature the higher minimum takes over. The temper- 
ature dependence of the magnetization in this region is 



(a) A = 



T/To 




M/Mo 



(b) A = -0.14J 




FIG. 4: The magnetization m in units of the T = sat- 
uration magnetization mo as a function of temperature at 
A = -0.14J. T = 17.43 K. The second order transition 
at T ~ 1.028 To, followed by the first order transition at 
T « 0.89 T . 
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FIG. 5: The magnetization m in units of the T = saturation 
magnetization mo in the paramagnetic phase at T = 12.55 K 
(T c = 12.2 K), as a function of magnetic field B, calculated 
from minimization of free energy for parameters t\ = 1.16 eV, 
t2 — 0.68 eV, J — 0.16 eV, appropriate to EuB(,. 



(c) A = -0.3J 

FIG. 3: Free Energy as a function of m at different tempera- 
ture, (a) A = where only first order phase transition occurs. 
From top to bottom: T > T cl , T = T cl , and T < T cl . (b) 
A = —0.14 J where both first and second phase transition 
occur. From top to bottom: T = T c2 , T c2 > T > T cl , and 
T < T c i. (c) A = —0.3 J her only second order phase tran- 
sition occurs. From top to bottom: T > T C 2, T = T C 2, and 
T < T c2 . 



shown in Fig|0J. Finally, panel (c) shows the free energy 
in the region in which there is no first order transition at 
all. We note that pressure, by increasing hybridization, is 
expected to broaden bands, thereby makes A more neg- 
ative. The predicted change in order of transition may 
therefore be accessible in pressure experiments. 

The presence, in the wide regions of the phase dia- 
gram, of a metastable high m free energy minimum, sug- 
gests that the system may undergo a metamagnetic phase 
transition. As shown in FigJ^J, this is indeed the case. 
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IV. APPLICATION TO EuB 6 



In this section, the model is applied to EuBq. 
The literature presents two conflicting interpretations 
of this compound: the semimetal interpretation fa- 
vored by quantum oscillations and bandstructure 
calculations [6j, and the large band gap interpretation, 
implied by photoemmsion|7j. These imply quite differ- 
ent parameters; we consider them seperately, beginning 
with the semimetal case. 

There are four parameters in this model - conduc- 
tion and valence bandwidths ti and t 2 , coupling strength 
J, and the bandgap A. t\ and t 2 are fixed by the 
effective masses of both bands given in |(J, and we 
found t\ = 1.16 el/, t 2 = 0.68 eV. Two sets of J 
and A are chosen - J — 0.14 eV, A = —0.2 J and 
J = 0.16 eF,-0.04 J < A < 0.04 J. The former choice 
are the best-fit parameters for fixed critical temperature 
and plasma frequency jump. However, these parame- 
ters produce a specific heat in disagreement with exper- 
iment and are in the region where the system has both 
1st and 2nd order transitions whereas experiment ap- 
parently yields only one transition Q We cannot fit T c 
and Auj 2 simultaneously in the 1st order transition re- 
gion. J = 0.16 eV, -0.04 J < A < 0.04 J are chosen the 
match the critical temperature, but the calculated Auj 2 
and ACy are roughly 2.5 and 2 times larger than the 
experimental data. 

With the two sets of parameters, we compute critical 
temperature T c , jump in plasma frequency Aw^, jump in 
specific heat ACy, core spin magnetization m, and latent 
heat AQ. The results are summarized and compared 
with experiment in the following table. 



Quantity 


Expt 


(0.14,-0.2J) 


(0.16, ±0.04J) 


Tc (K) 


12-14 


12.2 


12.2 ± 1 


AuS 2 (IQ'cm-' 2 ) 


1.625 


1.7 


4.3 ±0.2 


AC V 
(J/K per mole) 


12 


69.63 
(10.47^80.1) 


21.2 ±2.4 


to / too 


* 


0.3 

(0.28-^0.58) 


0.75 ±0.5 


AQ (J /Mole) 


* 


30.92 


90.5 ±2 



where the parentheses in the table represents ( J(eV), A). 
For the parameters J = 0.14 eV and A = -0.2 J, T c2 = 
16.1K. 

The existence of a metastable m ^ state implies 
the existence of a metamagnetic phase transition. For 
example, FigJSl shows M(B) calculated at a T slightly 
greater than T c . We take T = Yl.hhK (T c = 12.2K). We 
see that the jump occurs at B = 1.12 tesla, the magnetic 
field where the magnetization jump occurs as h* , Fig© 
plots h* as a function of temperature. 

Photoemission experiments at T > T c detects a wide 
range band gap (A — leV), and a fermi level about 
0, 2eV above the conduction band minimum. The car- 
riers in this band must arise from defects which dope the 
system; we model the defects as a flat band (7712 — > 00, 



in 

0) 




FIG. 6: h*(T), see the text for the definition of h* . The 
matemagnetic phase transition ends at T — 1.28 T c \ for J = 
0.16eV~, A = 0. The inset is for parameters J = 0.2eV", 
A = — 0.2J where the metamagnetic transition ends at T = 
1.012k . 
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FIG. 7: M(T) at photoemission parameters. T C 2 = 12.5K. 
The inset shows the square of plasma frequency (AtUp) (in the 
unit 10 7 cm -2 ) as a function of temperature. 



t 2 = 0). The measured conduction band dispersion 
implies t\ = 1.54eT/; the fermi level position yields 
A = -0.2eF; J 22 = 0, and by fitting T c we obtain 
J 11 = 0.14eV (quite close to the band theory value). 
These parameters place the system well into the second 
order region of the phase diagram in apparent contradic- 
tion to experiments. In this case, the jump in specific 
heat is 19.24 (J/Mole) while the plasma frequency and 
magnetization are all continuous as a function of temper- 
ature, as shown in fig'Q). 
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V. CONCLUSION 

We proposed a model to explain the combined 
metal/insulator, ferro/para magnetic phase transition of 
EuB 6 . The splitting of the Eu derived conduction band 
due to the coupling to core electrons at each Eu site is 
found to be the origin of the phase transition. From the 
technical point of view, the new feature of our results is 
the detailed examination of the dynamical mean field free 
energy. The free energy calculation implies the transition 
is first order. By choosing the physically reasonable pa- 
rameters, the available experimental data (T c , jump in 



plasma frequency, and jump in specific heat) are reason- 
ably well reproduced if the picture is favored by band 
theory and transported is adopted. The calculated jump 
in specific heat is around two times greater than the ex- 
periments. We predict a metamagnetic phase transition 
at fields up to 8 tesla, if the temperature is within a few 
degrees of T c \ . We predict that the modest pressures will 
change the order of phase transition from first to second. 
The parameters found by photocmission lead to a second 
order transition, a disagreement with other experiments. 

This work is supported by University of Maryland- 
Rut gers MRSEC. 
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